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Abstract. We have undertaken a series of hydrodynamic + N-body simulations in 
order to explore the properties of young stars. Our results suggest that the IMF 
may be sensitive to environment in its substellar region, with more brown dwarfs 
being formed where clusters are denser or more compact. We find that multiple 
stars are a natural outcome of collapsing turbulent flows, with a high incidence 
of A > 2 multiples. We find a positive correlation of multiplicity with primary 
mass but a companion frequency that decreases with age. Binary brown dwarfs 
are rarely formed, in conflict with observations. Brown dwarfs as companions are 
predominantly found orbiting binaries or triples at large separations. 
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1. Introduction and Motivation sible to perform calculations which both 

resolve the fragmentation and collapse of 
Most stars are known to be members of molecular clouds (accounting fully for the 
binary or even higher-order multiple sys- hydrodynamics) and which produce a sta- 
tems (Duquennoy & Mayor 1991). Thus, tistically significant number of steUar sys- 
any good star formation theory must be tems, thus opening the door to a direct 
a theory of (at least) binary star forma- comparison with observations of the IMF 
tion. Currently we can hope to do more and multiplicity properties of young stars 
than produce multiple stars by imposing (Delgado-Donate et al. 2004a,b). In this pa- 
some multi-armed instability on a coUaps- per we present the results from the first hy- 
ing core. Turbulent initial conditions, for drodynamic calculations to produce a sta- 
example, allow star formation to be trig- tistically significant number of stable mul- 
gered in a less predictable way (e.g. Bate et tiple systems in the separation range 1 - 
al. 2003). In addition, it has become com- IQOO AU. We will concentrate our attention 
putationally affordable to study the statis- mostly on those aspects pertaining to the 
tics of star pairing beyond simple N-body multiplicity properties of stars and brown 
integration (Delgado-Donate et al. 2003). dwarfs. 
These two steps forward have made it pos- 
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Fig. 1. Evolution of a typical a = — 3 cloud 
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Fig. 2. Mass Functions. Left, a = —3; right a = —5. The solid line includes all objects 
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2. IMumerical Scheme and Initial 
Conditions 

We have performed 10 calculations of small 
fragmenting clouds, using the SPH tech- 
nique to solve the fluid equations. Our 
version of SPH uses individual smoothing 
lengths and timesteps. Sink particles re- 
place bound blobs after a critical density is 
reached (Bate et al. 1995). We apply stan- 
dard viscosity with a. = 1 and /? = 2, 
and a binary tree to find nearest neigh- 
bours and calculate self-gravity. The opac- 
ity limit for fragmentation (see e.g. Low 

6 Lynden-Bell 1976) is modeled using a 
barotropic equation of state p cx p''', so 
that the gas is isothermal at low densi- 
ties (< 10~^^ g cm~^) and polytropic with 

7 = 5/3 at higher densities. 

Each cloud is initially spherical, has ra- 
dius of « 10^ AU, 5 Mq and density and 
temperature of « 10^^* g cm"'' and 10 K 
respectively. These values for p and T im- 
ply an initial Jeans mass of « 0.5 Mq, typ- 
ical of molecular clouds. We follow Bate & 
Burkcrt (1997) and use at least 100 SPH 
particles to resolve the minimum Jeans 
mass that can occur in the calculation (a 
few Mj), thus resulting in a total of 3.5 x 
10^ SPH particles. 

We impose a random 'turbulent' veloc- 
ity field on each calculation, defined by a 
power-law spectrum. The power-law expo- 
nent a is set to —5 in 5 of the simula- 
tions and to —3 in the other half (the for- 
mer index corresponds to shifting the bal- 
ance to having even more power in large 
scales than in the a = — 3 case). These val- 
ues of a bracket the observed uncertain- 
ties in Larson's velocity dispersion- size re- 
lationship. The velocity field is normalised 
so that there is cquipartition of kinetic and 
gravitational energy initially. The velocity 
field is allowed to decay freely. We are im- 
posing a parameterised initial velocity field 
which approximately reproduces observed 
bulk motions in molecular clouds (often 
described as 'turbulent' motions) but this 
term ('turbulence') should not be taken to 
imply that we are modeling what a fluid dy- 



namicist would recognise as fully developed 
turbulence. 

The gravitational force between sink 
particles is smoothed at short distances. 
Therefore, binaries with semi-major axis 
< 1 AU cannot form. A more detailed de- 
scription of the code and the initial condi- 
tions can be found in Delgado-Donate et al. 
(2004a). 

The calculations are run until star for- 
mation no longer occurs (see Figure 1 for 
snapshots of an a = —3 calculation). This 
translates into « 0.5 Myr of time or 60% ef- 
ficiency in terms of the amount of gas con- 
verted into stars. At this point the remain- 
ing gas is removed and the stellar system 
is evolved as a pure N-body system, using 
NB0DY2 (Aarseth 1999). After 10 Myr we 
find that 95% of the multiples have decayed 
into stable configurations (using the crite- 
rion by Eggleton & Kiseleva 1995), and we 
stop the integration. The calculations have 
been performed using the United Kingdom 
Astrophysical Fluids FacUity (UKAFF). 

3. Results: IMF 

Overall, the calculations produce 145 stars 
and brown dwarfs after 0.5 Myr, an average 
of « 15 objects per calculation (a number 
comparable to the initial number of Jeans 
masses). For the vast majority of objects, 
the accretion rates are very low and so one 
expects the objects masses to be represen- 
tative of their final masses. Overall, 40% 
of the objects are brown dwarfs and 60% 
stars. But the percentage changes depend- 
ing on the initial value of the turbulent 
spectrum slope. A higher fraction, nearly 
50% of the objects in the a = — 3 calcu- 
lations are brown dwarfs for » 30% in the 
a = —5 case. 

The MFs derived from the two sets of 
calculations can be seen in Figure 2. A 
KS test confirms that the MFs are dif- 
ferent at the substellar regime, at a 2a 
level, while the stellar MFs are indistin- 
guishable. Notice that the upper-end of the 
MFs cannot be directly compared with a 
Salpeter IMF as the MFs we present pro- 




ceed from calculations of clouds of the same 
mass, whereas the high-mass end of the 
observed IMF is probably a reflection of 
the cloud mass function, which is likely a 
power-law, non-flat distribution. It is natu- 
ral that the a = — 3 calculations produce a 
higher nmuber of brown dwarfs. These cal- 
culations have less kinetic energy in large 
scales relative to the a = —5 simulations 
and thus the action tends to concentrate in 
dense, compact clusters near the cloud cen- 
tre. This leads to much more frequent dy- 
namical interactions and consequently ejec- 
tion of low-mass objects, thus enhancing 
the fraction of brown dwarfs. Other initial 
conditions (e.g. a denser cloud, see Bate 
this volume) can also lead to a disparity in 
the fraction of brown dwarfs formed. Thus, 
the substellar IMF appears as the region 
of the IMF most sensitive to environment. 
Observational hints to this conclusion can 
be found in the literature, e.g. Briceno et 
al. (2002) find a paucity of brown dwarfs in 
Taurus relative to denser SFRs like Orion 
(Slesnick et al. 2004). 

4. Results: Multiplicity Properties 

Our simulations produce a wealth of mul- 
tiple systems. The multiplicity fraction at 
0.5 Myr after the initiation of star for- 
mation is close to 100%. The systems 
can adopt a variety of configurations, like 



binaries orbiting binaries or triples (see 
Figure 3). It is apparent that multiple star 
formation is a major channel for star for- 
mation in turbulent flows. 

4.1. Multiplicity as function of primary 

mass 

The companion frequency decreases during 
the first few Myr of N-body evolution, as 
many of the multiples are unstable. This 
internal decay affects mostly low-mass out- 
liers, which are released in vast amounts 
to the field. We expect that in a real clus- 
ter the multiplicity would drop even further 
as star forming cores do not form in isola- 
tion but close to one another. Some of our 
binaries orbiting binaries might not have 
survived in a more realistic environment. 
The predicted decrease in the multiplicity 
frequency has been quantitatively observed 
by Duchene et al. (2004). 

The properties of our multiples and 
the dependence of the binary fraction on 
primary mass is best illustrated by a di- 
rect comparison with the infrared colour- 
magnitude diagram of the 600 Myr old 
Praesepe cluster (Figure 4). This compar- 
ison summarises the conclusions obtained 
from a wider cross-check with observa- 
tions. The cluster was observed by Hodgkin 
et al. (1999) and our masses were con- 
verted to magnitudes using the tracks by 



E. Delgado and C. Clarke: VLMS Multiplicity Properties 



27 



10 



14 



16 



1.0 Mo 

• □ 

t 

m 

• 


Praesepe data • 
Model: multiples □ 


• 


Model: singles 

0.8 Mo 


. • • • 

#, CP □: 

• □ , 0.6 Mo 




" □ 

• 

□ 


error bar: 


g 0.4 Mo 
CP 


600 Myr 


• 

»m \ 

°n • 0.2 Mo 

«• • . 


Res. 2.5" (200 AU) 

o» . 

• 



1.0 1.5 2.0 2.5 3.0 

I-K 



Fig. 4. Colour-magnitude diagram (/ vs / — K) for the Praesepe clusters and superim- 
posed our results. Symbols as in legend. Resolution at 200 AU. 



Baraffe et al. (1998). Binaries with less than 
200 AU separation (the spatial resolution of 
Hodgkin et al.'s measurements) are consid- 
ered one object. 

Two features from Figure 4 are worth 
noting: first, the simulated cluster shows a 
binary sequence whose width is comparable 
to that of the Praesepe, except for systems 
redder than I — K = 2.5. This seems to 
suggest that the formation of a significant 
minority of triples, quadruples, etc. may in- 
deed be common in real clusters. Second, 
although our binary fraction for G stars is 
in agreement with observations, our models 
fail to produce as many low-mass binaries 
as observed. For example, at least a binary 



fraction of 15% is seen among brown dwarfs 
(e.g. Bouy et al. 2003). 

4.2. Where do we find brown dwarfs 

During the first few xlO'^ yr most brown 
dwarfs are locked in multiple systems, of- 
ten orbiting a binary or triples in eccentric 
orbits at large separations. Most of these 
systems are unstable and decay in a few 
Myr, releasing individual brown dwarfs to 
the field. Only a few substellar objects sur- 
vive bound to stars. Of these, the major- 
ity orbit a binary or triple at distances 
greater than 100 AU. One case out of 4 con- 
sists of a brown dwarf orbiting an M star 
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at 10 AU. Our results arc in agreement 
with the observed brown dwarf desert at 
very small separations (see e.g. Forveille, 
this volume). However, more than a dozen 
substellar objects companion to stars at 
wide separations arc known (Gizis ct al. 
2001). According to our results, we would 
expect that a large fraction of the primaries 
in these wide systems should turn out, in 
closer examination, to be > 2 multiples. 



5. Conclusions 

We have undertaken the first hydrodynam- 
ical + N-body simulations of multiple star 

formation to produce a statistically signifi- 
cant number of stable hierarchical multiple 
systems, with components separations in 
the range 1 — 1000 AU. We have shown 
that a high multiplicity fraction is typical 
of the very early stages, a few x 10''' yr after 
star formation begins, with many different 
possible multiple configurations. At later 
stages, a few Myr, most systems have 
decayed, ejecting brown dwarfs to the field 
and decreasing the companion frequency. 
Both the high initial multiplicity and its 
dependence with age seem to be in accord 
with resent observations. In addition we 
have probed different power spectra for 
the initial random velocity field and found 
that a larger fraction of brown dwarfs 
is produced when the initial conditions 
favour a more compact, dense distribution 
of stars. Our findings, taken together with 
others found in the literature, seem to 
suggest that the siibstellar regime is where 
it is most likely that the universality of the 
IMF might break down. 

We find a positive dependence of mul- 
tiplicity with primary mass, with few low- 
mass stars being primaries. The paucity of 
brown dwarf binaries in our simulations in- 
dicate that the models need finer tuning. 
Brown dwarfs are found, however, orbiting 
binaries or triples at large distances, and 
thus we suggest that a good test of our 
models is to look into the primaries of wide 



brown dwarf companions in search of mul- 
tiplicity. 
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